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ABSTRACT 

We report on the optical spectroscopic follow up observations of the candidate counterparts 
to four INTEGRAL sources: IGR J04069+5042, IGR J06552-1146, IGR J21 188+4901 and IGR 
J22014+6034. The candidate counterparts were determined with Chandra, and the optical obser- 
vations were performed with 1.5-m RTT-150 telescope (TUBITAK National Observatory. Antalya, 
Turkey) and 2.4-m Hiltncr Telescope (MDM Observatory, Kitt Peak, Arizona). Our spectroscopic 
results show that one of the two candidates of IGR J04069+5042 and the one observed for IGR 
J06552— 1146 could be active late- type stars in RS CVn systems. However, according to the likeli- 
hood analysis based on Chandra and INTEGRAL, two optically weaker sources in the INTEGRAL 
error circle of IGR J06552— 1146 have higher probabilities to be the actual counterpart. The candi- 
date counterparts of IGR J21188+4901 are classified as an active M-type star and a late-type star. 
Among the optical spectra of four candidates of IGR J22014+6034, two show Ho emission lines, one 
is a late-type star and the other is a M type. The likelihood analysis favors a candidate with no 
distinguishing features in the optical spectrum. Two of the candidates classified as M type dwarfs 
are similar to some IGR candidates claimed to be symbiotic stars. However, some of the prominent 
features of symbiotic systems are missing in our spectra, and their NIR colors are not consistent with 
those expected for giants. We consider the IR colors of all IGR candidates claimed to be symbiotic 
systems and find that low resolution optical spectrum may not be enough for conclusive identification. 

Subject headings: Stars: late-type - Stars: activity - Stars: flare - binaries: symbiotics - X- 
rays: individuals - X-ray:binaries 



1. INTRODUCTION 

The International Gamma-Ray Astrophysics Labora- 
tory (INTEGRAL; Winkler et al. 2003) has discovered 
hundreds of new hard X-ray sources (so called "IGR" 
sources) at energies above 20 keV since its launch on Oc- 
tober 2002. From over 700 sources that are listed in the 
4 th Bird Catalog (Bird et al. 2010), more than 400 of 
them are exclusively IGR sources 9 . INTEGRAL is very 
successful in finding new hard X-ray sources thanks to its 
wide field of view, observing strategy and most impor- 
tantly unique imaging capability of the ISGRI detector 
(Lebrun et al. 2003) with a localization accuracy of a 
few arcmin. However, such localization accuracy is often 
not enough to identify the correct optical counterparts of 
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these sources. Without multiwavclength information, it 
may not be possible to understand the physical origin of 
the X-ray emission and type of the emitting system. To 
identify a unique optical or infrared counterpart, arcsec- 
ond accuracy localization of the X-ray source is required. 

Our group observes the positions of IGR sources with 
soft X-ray telescopes such as Chandra and Swift, not only 
for identifying the correct counterparts, but also for pro- 
ducing 0.3—10 keV spectra that can be used to measure 
column densities and continuum shapes (Tomsick et al. 
2006, 2009, 2012; Rodriguez et al. 2008, 2009, 2010). 

From such spectroscopic and imaging efforts, it was 
possible to identify the nature of a large fraction of 
these IGR sources (Chaty et al. 2008; Butler et al. 
2009; Masetti et al. 2007, 2008, 2012, and references 
therein). Among these sources >250 are AGN. The 
Galactic sources include High-Mass X-ray Binaries, Low- 
Mass X-ray Binaries and Cataclysmic Variables (CVs), 
isolated neutron star (NS) systems as well as supernova 
remnants and pulsar wind nebulae. Yet, the nature of 
one third of IGR sources has not been identified. The 
identifications include a low number of unexpected hard 
X-ray sources like RS Canes Venatici (RS CVn), and 
symbiotic binary systems (Bird et al. 2010; Masetti et 
al. 2012; Rodriguez et al. 2010), and thus searches for 
more such objects are warranted. 

So far, there have been four IGR sources tentatively 
identified as RS CVn binaries (Rodriguez et al. 2010; 
Masetti et al. 2012; Krivonos et al. 2010), but no firm 
association with RS CVn systems and IGR sources exists. 
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In general, RS CVn systems are the members of chromo- 
sphcrically active stars which include flare stars (dMe 
stars or UV Ceti stars) and T Tauri stars as well. These 
systems are detached binaries having emission lines of 
Call H & K and Ha lines in their optical spectra. The 
components of an RS CVn system can be G-M type stars 
within the luminosity classes of II- V (Kogure & Leung 
2007). This type of system is known to have strong X- 
ray and radio emission (Barbier et al. 2005; Osten et al. 

2007) . They show solar like magnetic activities and flares 
(Haisch et al. 1990). 

Symbiotic systems are long period interacting binaries 
consisting of a red giant star (G-M) and a hot compan- 
ion surrounded by an ionized nebula (Friedjung & Viotti 
1982; Kenyon 1986, 1990; Kenyon & Webbink 1984). The 
hot component is most often a white dwarf (WD) or a 
disk-accreting main sequence star. On the other hand 
there are also systems having a NS as the compact com- 
panion (Mikolajewska et al. 1997). In general, the optical 
spectra of symbiotic stars have photospheric absorption 
features and molecular bands of the cool giant and emis- 
sion lines of the HI Balmer series, Hel, Hell, [OIII] and 
[Nelll] which are produced in the surrounding nebula of 
the hot component by its intense UV radiation (Kenyon 
&: Fernandez- Castro 1987; Cieslinski et al. 1994). In ad- 
dition, these systems can show hard, soft and supersoft 
X-ray emissions (Kennea et al. 2009; Muerset et al. 1997). 

So far, eleven IGR sources have been suggested to 
be symbiotic systems according to their optical/infrared 
features and X-ray properties. Among those identi- 
fied, the symbiotic nature of IGR J10109-5746 (CD-57 
3057; Masetti et al. 2006; Kennea et al. 2009), IGR 
J12349-6434 (RT Cru; Masetti et al. 2005; Luna 
& Sokoloski 2007; Kennea et al. 2009) and IGR 
J16194-2810 (Masetti et al. 2007; Ratti et al. 2010) are 
well determined through multiwavelength analyses. A 
nearby system (1.56 kpc, parallax measurement), IGR 
J15293— 5609, has been classified as a symbiotic binary 
based on the well determined surface temperature and 
radius of the companion star (Tomsick et al. 2012). On 
the other hand, the identifications of IGR J16358— 4726 
(Nespoli et al. 2010), and IGR J17497-2821 (Paizis et 
al. 2007, 2009) as symbiotic systems are questionable. 
Chaty et al. (2008) suggest that IGR J16358-4726 could 
be a high mass X-ray binary. Likewise the classifica- 
tion of IGR J17497— 2821 is also controversial due to 
the suggestion of a black hole primary (Walter et al. 
2007; Paizis et al. 2007) would make the system very un- 
usual among the known symbiotics. IGR J16393— 4643 
was also claimed to be a symbiotic system based on the 
.ftf-band spectrum of the candidate counterpart 2MASS 
J16390535-4642137 (Nespoli et al. 2010) which was pre- 
viously classified as a BIV-V type star using its opti- 
cal/IR spectral energy distribution (Chaty et al. 2008). 
Later, it was revealed that the refined X-ray position of 
IGR J16393— 4643 was incompatible with this controver- 
sial 2MASS source (Bodaghee et al. 2012). Candidate 
counterparts to IGR J11098-6457, IGR J17197-3010 
and 1RXS J174607.8-213333 (Masetti et al. 2009, 2012, 

2008) are also tentatively claimed to be late M type gi- 
ants based on comparison of their optical spectra with 
the known symbiotic systems, but it was later shown that 
the tentative optical counterpart identifications of IGR 
Jl 1098-6457 and IGR J17197-3010 are wrong (Tomsick 



et al. 2009; Luna et al. 2012). Finally, IGR J16293-4603 
was suggested to be a candidate symbiotic system by 
Ratti et al. (2010) based solely on its optical/NIR pho- 
tometry. 

In this work, we investigated the optical candidate 
counterparts of four IGR sources in the 4 th Bird Cat- 
alog (Bird ct al. 2010) with unknown X-ray classi- 
fication: IGR J04069+5042, IGR J06552-1146, IGR 
J21188+4901 and IGR J22014+6034. The fields of these 
sources were observed with Chandra, and possible opti- 
cal/infrared counterparts were identified (Tomsick et al. 
2012, hereafter T12). All four sources have low Galac- 
tic latitudes (see Table 1), and high extinction limits the 
possible counterparts mostly to relatively nearby stellar 
objects. We performed optical spectroscopy of bright 
candidates in the USNO and 2MASS catalogs, and also 
obtained images that were located near the brightest 
Chandra sources. This work is follow up of T12 that uti- 
lized Chandra to find X-ray counterparts of IGR sources 
for further investigations. 

The paper is organized as follows. In Section 2, we 
summarize the observations and data reduction in the 
optical, in Section 3 we identify the types of optical coun- 
terpart candidates of the IGR sources, and in the final 
section, we discuss the possible counterparts to each IGR 
source, and compare the M type candidates we found 
with earlier tentative symbiotic identifications. 

2. OBSERVATIONS AND DATA REDUCTION 

The INTEGRAL coordinates and the fluxes of the four 
sources we investigated are given in Table 1. The Chan- 
dra localizations (See Table 2) of the candidate counter- 
parts allowed us to select sources to be observed at the 
TUBITAK Turkish National Observatory 10 (TUG) and 
at the MDM Observatory. A detailed list of all detected 
possible Chandra counterparts and their soft X-ray prop- 
erties are reported in T12. All of the X-ray sources de- 
tected by Chandra in the field of each IGR source are 
shown in Figure 1. We also searched the online cata- 
logs of 2MASS, USNO-B1.0 and USNO-A2.0 to find the 
IR and optical counterparts associated with the Chandra 
positions (See Table 3). 

The spectroscopic observations of candidate counter- 
parts for IGR J04069+5042 and IGR J06552-1146 were 
done at the MDM Observatory, while IGR J21188+4901 
and IGR J22014+6034 were observed at TUG. 

2.1. TUG Observations 

The medium-resolution spectra of the candidate coun- 
terparts to the sources IGR J21 188+4901 and IGR 
J22014+6034 were obtained with the TUBITAK Faint 
Object Spectrometer and Camera (TFOSC) which 
is mounted on the Russian- Turkish 1.5 m Telescope 
(RTT150) located at TUG, Antalya, Turkey. The cam- 
era is equipped with a 2048 x 2048, 15 /mi pixel Fairchild 
447BI CCD. It has a FOV of 13' x 13' with a pixel scale 
of 0'.'39 pixel -1 . Grism #15, having an average disper- 
sion of 3 A pixel -1 , and slit 67 /im (1'.'24) were used 
providing a 3300-9000 A wavelength band. 

All spectroscopic data, acquired on 2011, August 27 
under good weather conditions, were reduced using the 

10 http : //www . tug . tubit ak . gov . tr 



Investigating the Optical Counterpart Candidates of INTEGRAL Sources 



3 



Long-Slit package of MIDAS 11 to obtain ID spectra of 
target sources. The spectra were corrected for bias, flat- 
fielded and cleaned from the cosmic-ray hits with stan- 
dard MIDAS routines. Wavelength calibration of the 
spectra were done by Neon lamp spectrum. The helio- 
centric correction was applied to each spectrum after the 
extraction was done. The spectrophotometric standard 
star BD +33°2642 (Oke 1990) was observed during the 
observing night to get the flux-calibrated spectra of the 
candidate counterparts. 

The equivalent widths (EW) of prominent absorption 
and emission lines are measured by the ALICE subrou- 
tine of MIDAS. 

2.2. MDM Observations 

We used the OSMOS (Ohio State Multi-Object Spec- 
trograph) with the MDM4K CCD on the 2.4 m Hiltncr 
Telescope of the MDM Observatory on 2011 January 11 
and 13. The Volume Phase Holographic (VPH) trans- 
mission grating we utilized has a l'/2 wide slit with a 
wavelength coverage of 3900-6800 A at 3.5 A resolution. 
Conditions were clear during the observing runs. Spec- 
tra of two Chandra localized X-ray sources in the IN- 
TEGRAL error circle of ICR J04069+5042, and of one 
X-ray source for IGR J06552— 1146 were obtained. The 
spectral reduction was performed using standard IRAF 12 
routines. The flux calibration was done by Oke & Gunn 
(1983) standard stars. We also obtained images in the 
i?-band using the same CCD for sources that are not in 
the standard catalogs. 

2.3. Chandra Observations 

We used the Chandra observations not only for de- 
termining optical counterpart candidate positions, but 
also obtaining soft X-ray properties which could aid 
in the identification process. We fitted the Chandra 
0.3-10 keV spectra with an absorbed power-law model 
to determine the absorbed 0.3-2 keV and 2-10 keV 
flux for each source. We defined the hardness ratio as 

(-Fa-lOkeV — ^0.3-2kev)/(-F2-10keV + ^0.3-2kev)- F° r de- 
tails of spectral fitting, background subtraction and flux 
calculation, see T12. For soft X-ray properties relevant 
for this work, see Table 2. 

3. RESULTS 

3.1. IGR J04069+ 5042 

We performed optical spectroscopy of two possible 
candidates localized by Chandra within the INTEGRAL 
error circle of 4.'2. 

2MASS J04064392+5044469 (#1): The optical spec- 
trum includes prominent absorptions of Nal D resonance 
lines, Mgl b and Cal triplet lines (See Fig. 2). Although 
distinguishing features in the classification region are 
mostly at the noise level, the spectrum is very similar to 
the K type main sequence star HD 283916 (SAO 76803) 
in Jacoby et al. (1984) and to the sources within the 
same class in Ross et al. (1996). Indeed the presence 
of MgH bands at around A4780 and A5200, seen only 
in the spectra of main sequence stars, supports the 

11 http : //www . eso . org/pro j ects/esomidas/ 

12 http://iraf.noao.edu/ 



luminosity classification. On the other hand a weak but 
noticablc Ha emission line (EW = 1.07 A, sec Table 4) 
which could be due to the chromospheric activities, 
stellar winds and/or interactions with the environment 
of the secondary (or primary) indicates an active star. 
Therefore we classify this source as a K5-7 Ve type star 
based on its spectral lines and EW measurements. 

2MASS J04064872+5039316 (#2): The spectrum 
is dominated by neutral Hydrogen lines (H<5, H7, H/3 
and Ha) and several metallic lines of Fel, SrII, Mgl and 
MgH (See Fig. 3). The Cal triplet lines that are used 
as luminosity indicators are very weak (See Table. 4) as 
seen in the spectra of F type stars (Montes & Martin 
1998; Cayrel et al. 1996). In addition the G band of 
CH at A4300 that becomes visible for the spectral types 
later than F4 and increases in strength with decreasing 
temperature (Gray & Corbally 2009) is one of strongest 
features in the blue region (4000-5000 A) and allows us 
to exclude the earlier spectral classes. Comparing the 
line width ratio of A4077 (SrII) to A4226 (Cal) that is 
a sensitive criteria for luminosity (Morgan et al. 1943), 
and the EW measurements of Balmer lines (See Table 4) 
to the values in Jaschek & Jaschek (1995) and in Eaton 
(1995) we suggest that this source could be an F5-8 
type star within the luminosity class of III-V. 

3.2. IGR J06552-1146 

Optical spectroscopy was only obtained for candidate 
#1 (see Fig. 1) which is ~<y/8 away from its Chandra 
counterpart. Candidate #2 is 2'.'8 away from the 
nearest source in 2MASS and USNO catalogs. As it 
can also be seen in the i?-band image of this source 
(Fig. 4), the Chandra source cannot be clearly identified 
with an optical counterpart. Likewise, candidate #3 
was not listed in Table 1 since its coordinates are 
not associated with any optical or IR counterparts in 
the catalogs. However, a weak source in the Chandra 
error circle is visible in the i?-band image of the field 
(Fig. 4). This source has an R magnitude of 22.17 ± 0.14. 

2MASS J06545833-1149119 (#1): The flux-calibrated 
spectrum of candidate #1 shows Fel, Mgl and Cal 
absorption lines (See Fig. 5) and a weak (EW = 0.45 
A) Ha emission line (see Table 4). The characteristics 
of the spectrum are very similar to those of candidate 
#1 to IGR J04069+5042 in terms of the existence of 
the spectral lines and the molecular MgH band. To give 
a spectral range, we searched the spectral atlases and 
libraries having a similar spectrum to that of 2MASS 
J06545833— 1149119. By comparing the spectra of the 
sources listed in Jacoby et al. (1984), Ross et al. (1996) 
and in Gray & Corbally (2009) to candidate #1 we 
suggest that it would be a main sequence K5-8 type 
emission-line star. 

3.3. IGR J21188+4901 

This source is tentatively classified as a transient in 
Bird et al. (2010) due to its highly variable X-ray flux in 
20-40 keV energy band. Two candidates which are quite 
faint in optical band have been observed (Fig. 1) to 
search for the optical counterpart to IGR J21188+4901. 
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2MASS J21182288+4906259 (#1): The optical spec- 
trum of this source contains both broad molecular TiO 
bands with a red continuum and an Ha emission line 
with EW = 1.9 A. Although EW values of H 7 and H/3 
lines could not be determined due to their low signal 
to noise, they are apparent in emission. These spectral 
features are indicative of a late-type star spectrum with 
chromospheric activity. Cal and Call triplet lines, used 
as luminosity indicator for M type stars (Jaschek & 
Jaschek 1995), Mgl lines and telluric absorption bands 
are the other identified features in the spectrum (See 
Fig. 6). We point out that the spectral apperance of 
this candidate is very similar to that of the counterparts 
of 1RXS J174607.8-213333 and ICR J11098-6457 
which arc identified as an M2 4 (Masetti et al. 2008) 
and an M2-type (Masetti et al. 2009) symbiotic giants 
respectively. However this source does not have the 
clear spectral features of a red giant star. In addition 
2MASS magnitudes are too faint to be a nearby giant. 
Therefore we searched for the late-type stars spectra 
in the atlases of Gunn & Stryker (1983), Jacoby et al. 
(1984) and Silva & Cornell (1992) to compare with our 
flux-calibrated spectrum and we determine that the 
spectral class of this candidate could be a M3-5 type 
active main sequence star. 

2MASS J21183906+4858049 (#2): This candidate 
is the only source having a Chandra counterpart within 
the INTEGRAL error circle. The spectrum shown in 
Fig. 7 has typical features of a late-K/mid-M type star 
(Gray & Corbally 2009) without any emission lines. 
In addition the molecular CaOH band that becomes 
visible for dwarfs having spectral class later than K7 
and MgH bands are prominent features in the spectrum. 
The weakness and/or absence of the Balmer lines (See 
Table 5) support the late-type spectral classification. 
Therefore, considering characteristics of the spectrum 
we suggest a K7-M4 type main sequence classification 
for this source. 

3.4. IGR J22014+6034 

We obtained optical spectroscopy of four stars that are 
detected with Chandra in 2MASS and USNO catalogs. 
Two of these candidate optical/IR counterparts lie 
within the 5'.4-radius INTEGRAL error circle as seen in 
Figure 1. 

2MASS J22002116+6033420 (#1): The flux-calibrated 
spectrum is very similar to candidate #2 of IGR 
J21188+4901 that is classified as a K7-M4 type main 
sequence star. It is dominated by the molecular MgH, 
CaOH and CH-G bands in addition to Mgl b triplet 
lines (See Fig. 8). The emission lines of Ha and Call 
H & K at A3968, A3934 (See Table 5), indicators of 
chromospheric activity, are also detected. 

2MASS J22010934+6034088 (#2): This IR source 
is the optically brightest candidate to the X-ray position 
of IGR J22014+6034. However, the optical spectrum 
contains no significant features related to an X-ray 
source. The most prominent absorptions in the blue 
region of the spectrum are of Call H & K and Balmer 
series lines. Metallic lines of Fel and Mgl b triplet are 
also visible. The flux-calibrated spectrum of 2MASS 



,122010934+6034088 is similar to that of A-G type 
stars due to its decaying blue continuum towards the 
longer wavelengths (sec Fig. 9). On the other hand the 
presence of CH-G band, and the weakness of Balmer 
lines (See Table 5) while compared to A-type stars 
indicate that this source must be in F-G spectral range. 
By comparing EW values of Call triplet lines in near 
IR region to those of late-type stars (Jaschek & Jaschek 
1995; Zhou 1991), we classify this candidate as an F5-G5 
type III-V star. 

2MASS J22015416+6038094 (#3): This bright IR star is 
located at edge of the IGR error circle. It does not have 
an optical counterpart in USNO-A2.0 catalog. Although 
the flux-calibrated spectrum contains nearly most of the 
features detected in 2MASS J22010934+6034088, they 
arc different regarding the strength of the lines (See 
Table 5). The identification spectrum of the candidate 
#3 is shown in Figure 10. The Balmer lines of H<5 and 
H7 are extremely weak or absent, compared to Ha line 
(EW = -2.3 A) and to CH-G band. The blue region 
of the spectrum which is dominated by metallic lines 
of Mgl, MgH and Fel displays typical spectral features 
of a late-type star. We derive a spectral class of G 
III-IV based upon the prominent features for this source. 

2MASS J22020837+6030425 (#4): The coordinates of 
the source are coincident with the positions of USNO- 
B1.0 1505-0322924 and USNO-A.2 1500-08693376 
although they fall just outside the INTEGRAL error 
circle (6'06). The flux-calibrated 4000-9000 A spectrum 
shows typical features of a M type star (see Fig. 11) 
dominated by strong molecular bands of TiO (Gray & 
Corbally 2009). The emissions of Balmer series lines 
(H7, H/3, Ha) superposed with TiO bands and the 
presence of Cal, Call, Mgl b lines are in favor of the 
late-type spectral classification (See Table 5). Since 
determining exact spectral class of late-type stars is 
rather complicated due to the contaminated blue part 
of the spectrum, we used spectrophotometric atlases of 
Gunn & Stryker (1983), Jacoby et al. (1984) and the 
library of Silva & Cornell (1992) to secure the spectral 
and luminosity classification. Comparing our spectrum 
with the late-type stars in the atlases and libraries we 
classify 2MASS J22020837+6030425 as a main sequence 
M2-5 type chromosphcrically active star. 

3.5. X-ray results 

The four IGR sources we investigated in this work 
are reported as sources for which none of the Chandra 
sources in the field of view stands out as a candidate 
counterpart from a likelihood analysis just based on X- 
rays as described in T12. This likelihood analysis calcu- 
lates the spurious association probability as: 

P = 1 _ e -JV(>-F'2-10keV) T el earch Q\ 

where ^(>F 2 _ 10ke v) = 9.2 (F 2 _ lokcV /10- 13 )-°- 79 deg" 2 
(Sugizaki et al. 2001, T12), and i*2-iokcV is the absorbed 
2-10 keV flux in ergs cm -2 s _1 . 9 sea rch is the radius of 
the search region in units of degrees and equals to 90% 
confidence for the sources inside the INTEGRAL error 
circle (Sec Tabic 1). For sources outside the error circle, 
it indicates the angular separation from the center of 
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the INTEGRAL error circle. The spurious association 
probabilities of all sources are given in Table 2. 

For most of the candidates the spurious association 
probability is very high. The notable exceptions are can- 
didate #3 of ICR J06552-1146, candidate #2 of IGR 
J21188+4901, and candidate #3 of IGR J22014+6034. 
However the likelihood analysis does not take into ac- 
count the variability. The soft X-ray measurements were 
obtained at an arbitrary epoch, and INTEGRAL analysis 
shows that the maximum hard X-ray flux and the average 
hard X-ray flux are different for these sources indicating 
variability (See Table 1, T12, Bird et al. 2010). In fact, 
IGR J04069+5042 and IGR J06552-1146 have burstic- 
ity 13 impact grater than 1.1 and for IGR J21888+4901 
it is greater than 4, indicating a highly variable nature 
of the source in X-rays (Bird et al. 2010). Therefore like- 
lihood analysis based on X-ray data alone may fail in 
revealing the true counterpart if the X-ray observations 
take place when the source is weak. It is still possible 
that one of the detected Chandra sources is the correct 
counterpart. 

4. DISCUSSION 

Below, we discuss the candidate counterparts for each 
IGR source in terms of their optical spectra, Chandra 
soft X-ray properties and relative probabilities to find 
out the candidates that are more likely to be the IGR 
source. 

4.1. IGR J04069+ 5042 

According to the Equation 1, the spurious association 
probabilities of candidates #1 and #2 are >68%, and > 
99.99% respectively, and neither of them was detected 
in the 2-10 keV band with Chandra. Among these two 
sources, candidate #2 has typical features of a late F- 
type star. On the other hand, candidate #1 shows H Q 
emission, which makes it more likely to be the correct 
counterpart. 

4.2. IGR J06552-1146 

Candidate #1 is 4.'14 away from the center of the IN- 
TEGRAL error circle and the probability of spurious as- 
sociation is 88%. The spectrum shows evidence of an 
H Q emission line and it is similar to that of candidate 
#1 for IGR J04069+5042. The other sources with high 
ACIS count rates are too dim in the optical to get a 
spectrum. Their spurious association probabilities are 
33-98% and 14-29% for candidates #2 and #3 respec- 
tively. Candidate #3 in fact has positive hardness ratio, 
and the relative probability is quite close to the cut-off 
for a Chandra/ INTEGRAL association in T12. 

4.3. IGR J21888+4901 

Candidate #1 is outside the 90% confidence error circle 
of INTEGRAL with spurious association probability of 
>99.96%. It has been detected with Chandra in the 2-10 
keV band, and is an active M type star. Candidate #2 is 
the only soft X-ray source (detected in 0.3-2 keV energy 
band) in the INTEGRAL error circle with a probability 

13 Bursticity can be defined as the ratio of the maximum signifi- 
cance on any timescale, compared to the average significance (Bird 
et al. 2010). 



of >15% but it is not detected with Chandra in the 2-10 
keV band. Candidate #2 shows properties of a late-type 
(K7-M4) main sequence star. Based on probabilities, 
candidate jj=2 is a more likely counterpart than candidate 
#1- 

4.4. IGR J22014+6034 

The first candidate is the brightest Chandra source in 
the vicinity of the IGR source, however it is out of the 
90% confidence error circle of INTEGRAL. The spec- 
trum is typical of an active late-type star. The probabil- 
ity of spurious association is >95% and it is detected in 
the 2-10 keV band with Chandra, however the source is 
very soft. 

The optical spectra of candidates #2 and #3 show that 
they are late-type main sequence stars without emission 
lines. The probabilities of spurious associations for these 
sources are >88% and >36% respectively. 

The optical spectrum of candidate #4 shows features 
of an active M type main sequence star. It is located just 
outside the error circle of INTEGRAL. The probability 
is >99.98%, and it is not detected with Chandra in the 
2-10 keV band. 

4.5. IGR sources identified as symbiotic systems 

At an early phase of our analysis, we compared our 
optical spectra with those published in Masetti et al. 
(2008, 2009, 2012) and realized the similarities of candi- 
date #1 of IGR J21888+4901 and candidate #4 of IGR 
J22014+6034, the M type stars with broad TiO bands, 
to those identified as symbiotic stars. However, these 
sources do not show Hel, Hell and [OIII] in emission 
(Kenyon 1986), and their infrared colors are not con- 
sistent with the colors of symbiotics in the J — H vs 
H — Ks diagram as given by Phillips (2007) and Cor- 
radi et al. (2008). Our candidates are more probably 
active main sequence stars e.g., dMe stars (Pettersen & 
Hawley 1989) rather than being giants. To illustrate this 
we plotted the NIR color-color diagram of M type TUG 
candidates (filled red circles) along with colors of several 
other sources in Figure 12. The plot includes validated 
IGR symbiotics (IGR J10109-5746, IGR J12349-6434 
and IGR J16194-2810 with filled blue circles) and the 
all known symbiotics (small, empty, purple circles), the 
boundaries that limit stellar (S) and dusty (D) type sym- 
biotics (see Pereira & Roig 2009, for details of S and D 
type symbiotics) and the loci of main sequence stars and 
red giant stars. The colors of IGR sources, tentatively 
classified as symbiotic systems but later shown to have 
the incorrect counterpart candidates, are also placed in 
Figure 12. The initial classifications of the counterparts 
to IGR J11098-6457, IGR J17197-3010 (empty red cir- 
cles) and IGR J16393— 4643 (empty orange circle) were 
done with optical and/or IR spectral analyses. Since 
the refined X-ray positions of these IGR sources were 
not related to the suggested counterparts, the symbiotic 
identifications were ruled out. Note that the colors of 
incorrect counterparts can easily be associated with the 
main sequence stars like the TUG candidates. 

Two other candidates were also tentatively claimed 
to be symbiotic binaries, IGR J16358-4726, IGR 
J17497— 2821, but subsequent analyses made the sym- 
biotic interpretation questionable (see Introduction for 
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details). These sources are shown as "IGR Sym, ques- 
tionable" in Fig. 12 (filled brown circles), and have colors 
consistent with dusty symbiotic systems. The colors of 
the infrared counterpart of IGR J16293— 4603 (filled or- 
ange circle) are also consistent with being a symbiotic 
system as expected, since the initial identification was 
done using g, r, and i colors by Ratti et al. (2010). We 
also show IGR J15293-5609, classified symbiotic 
with a K type companion (blue star). This source falls 
slightly out of the region of validated symbiotics. On 
the other hand, symbiotics with K type companions are 
rare, and occupy lower J — H values (Pereira & Roig 
2009) just like IGR J15293-5609. These are called yel- 
low symbiotic systems which are indicated separately in 
Figure 12 (filled, dark yellow circles). We also note that 
reddening can move the sources along the diagonal lines 
shown in the figure, and the large uncertainty in Ay of 
IGR J15293-5609 (T12) could easily move this source to 
the confirmed region. Finally we show the colors of all 
other candidates that we identify as main sequence stars 
(with triangles) , and as expected most of them lie along 
the locus of main sequence stars. 

This analysis shows that near infrared color diagram 
provides a quick method to double check identifications 
of symbiotics with some caveats. First of all, many sym- 
biotic systems have normal NIR colors (see confirmed 
symbiotics, and the locus of red giant branch stars in 
Fig. 12). NIR colors may be enough to claim a giant 
identification, but it does not mean that the system is 
a symbiotic. Second, most Galactic IGR sources are in 
the plane with high extinction (including the sources we 
analyzed, all with Ay > 3) which could move main se- 
quence stars into the range of symbiotics. We could not 
include colors of 1RXS J174607.8-213333 (Masetti et al. 
2008) on this plot because it was not in the near infrared 
catalogs, but that identification also relied only on low 
resolution optical spectrum which requires further anal- 
ysis. As mentioned in Masetti et al. (2008), it could also 
be a spurious detection. 

5. CONCLUSIONS 
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TABLE 1 

The X-ray measurements of IGR sources given in Bird et al. (2010) 



IGR Name 


RA(J2000) 


DEC(J2000) 


Z a 


b h 


Error R,adius c 


FluX 2 0-40 keV 


Flux 40 _ 10 keV 


Peak Flux d 2o-40 


J04069+5042 


04 h 06 m 55?0 


+50°42.'l 


151.43 


-1.03 




1.6+0.3 


<1.0 


2.0+0.3 


J06552-1146 


06' 1 55 m 10?0 


-11°46.'2 


223.85 


-4.52 


its 


1.1+0.3 


<1.0 


2.6+1.0 


J21 188+4901 


21 h 18 m 48?0 


+49°01.'0 


91.27 


-0.33 


4.' 3 


<0.2 


<0.4 


4.5+1.3 


J22014+6034 


22 h 01 m 27 s 


+60°34.'0 


103.49 


+4.28 


5'A 


<0.2 B 


<0.4 





Note. — The fluxes are given in units of mCrab. A flux of 1 mCrab in the 20-40 keV energy range corresponds to 7.57x10 1 erg cm 
s" 1 while it is 9.42xl0~ 12 erg cm -2 s" 1 for 40-100 keV band. 
a The Galactic Longitude of the source in degrees. 
b The Galactic Latitude of the source in degrees. 
c The radius of the 90% confidence IGR circle. 

d The peak flux of the sources which show variability in the 20-40 keV energy band. For IGR J22014+6034 no variation detected in this 
range. 

e The source is detected in 17-30 keV band, Bird et al. (2010). 
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TABLE 2 

Chandra Localizations, Counts, Hardness Ratio, and The Probability of Spurious 

Associations 









IGR 


J04069+5042 








Candidate^ Chandra Name" 


RA(J2000) 


DEC(J2000) 


ACIS Counts 


Hardness b 


P (%) c 


^XOU J040643. 9+504446 


04' 1 


06 m 


43? 99 


+50°44' 


46'.' 2 


5.4 


-1.19+1.10 


>68 


2 CXOU J040648. 7+503931 




06 m 


48? 69 


+50° 39' 


3r72 


?,A 


-1.30+1.65 


>99.99 








IGR 


J06552 


-1146 








lPVnTT TOfl^/I^Q Q 11/1011 

JUD04DO.O — 114yil 


06 h 


54 m 


58?31 


-11°49' 


11'.' 3 


O . • J 


— U.l 0±U. / u 


>oo 


tAUU JUDOOZo.D — 1140U1 


06 h 


55 m 


23? 65 


-11°46' 


01"1 


1 Q A 


— U.DZ±U.4D 


QQ no 
oo-yo 


3 CXOU J065529. 5-114900 


06 h 


55 m 


29? 56 


-11°49' 


OCX'! 


35.5 


+0.02+0.22 


14-29 








IGR 


J21188+4901 








^XOU J211822. 9+490627 


21 h 


18 m 


22? 86 


+49° 06' 


26'.' 7 


7.1 


-0.78+0.82 


>99.96 


2 CXOU J211839. 1+485806 


21'' 


lg m 


39? 09 


+48° 58' 


05"5 


4.3 


+0.61+1.00 


>15 








IGR 


J22014 


+6034 








^XOU J220021. 1+603342 


22 h 


00 m 


21? 14 


+60° 33' 


41'.' 8 


45.8 


-40.91+0.25 


>95 


2 CXOU J220109. 3+603409 


90': 


01 m 


09?34 


+60° 34' 


08'.'7 


14.4 


-0.79+0.46 


>88 


3 CXOU J220154. 1+603809 


22': 


01 m 


54? 13 


+60° 38' 


09'.'5 


31.7 


-0.61+0.27 


>36 


4 CXOU J220208.4+603042 


22': 


02 m 


08? 35 


+60° 30' 


42'.' 4 


10.7 


-1.15+0.70 


>99.99 



a All Chandra information is reported by T12 

b Hardness ratio is defined as (C2— lofceV ~ Co. 3— 2kev) / (C2— WkeV + Cb.3-2fceV~) where C is the count 
rate in the given band. After background subtraction in Chandra, the hardness ratios may become 
negative (T12). 
c See Equation 1 and T12. 
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TABLE 3 

Observed Possible Optical/Ir Counterparts To The Sources 



Candld<1 ' e # Catalog/Source Name 3 


RA(J2000) 


IGR J04069+5042 

DEC(J2000) Distances 13 




Magnitudes 




1 2MASS J04064392+5044469 
USNO-A2.0 1500-08672098 
USNO-B1.0 1505-0322270 


04 ,l 06 m 43!92 
04 h 06 m 43?94 
04 h 06 m 43?92 


+50°44'46"95 
+50°44'47('39 
+50°44'47?32 


3f20/l'/00 
3?23/l'/28 
3.'21/1'.'29 


Bi= 


J=14.40 H=13.68 K s =13.45 
B=18.60 R=16.60 
=19.06 Ri ==16.32 B 2 ==18.49 R 2 ==16.71 I== 


=15.42 


2 2MASS J04064872+5039316 
USNO-A2.0 1350-04221566 
USNO-B1.0 1406-0110577 


04 h 06 m 48?73 
04 h 06 m 48?74 
04 h 06 m 48?73 


+50°39'31'.'63 
+50°39'31^'81 
+50°39'31'/76 


2.'76/0 / .'58 
2?75/0'.'78 
2. '76/0'.' 69 


Bi = 


J==11.44 H==11.13 K s ==11.06 
B=13.20 R=12.30 
=13.56 Ri==12.15 B 2 ==13.37 R 2 ==12.69 I== 


= 11.82 


IGR J06552-1146 


1 2MASS J06545833-1149119 
USNO-A2.0 0750-03064876 
USNO-B1.0 0781-0141882 


06 h 54 m 58!34 
06 h 54 m 58!35 
06 h 54 m 58!33 


-11°49'12^'00 
-11°49'11'.'80 
-11°49'11'.'89 


4.' 14/0'.' 79 
4.' 14/0'.' 74 
4.'14/0'.'64 


Bi = 


J==14.04 H==13.83 K s ==13.32 
B=16.60 R=15.40 
=16.95 Ri ==15.43 B 2 =17.48 R 2 =15.34 1= 


=14.60 


IGR J21188+4901 


1 2MASS J21182288+4906259 
USNO-A2.0 1350-13821833 
USNO-B1.0 1391-0394771 


21 h 18 m 22!89 
21 h 18 m 22?89 
21 h 18 m 22!88 


+49°06'25^'96 
+49°06'26^'39 
+49°06'26^'31 


6'.Sl/(f.'78 
6.'82/0'.'40 
6.'82/0'.'44 


Bi = 


J=13.01 H=12.45 K s =12.17 
B=18.20 R=16.30 
=18.33 Ri=16.01 B 2 =18.24 R 2 =16.06 1= 


= 14.51 


2 2MASS J21183906+4858049 
USNO-A2.0 1350-13829476 
USNO-B1.0 1389-0393707 


21 h 18 m 39?07 
21 h 18 m 39?00 
21 h 18 m 39?06 


+48°58'04('92 
+48°58W32 
+48°58W28 


3.'27/0 / .'61 
3.'26/0'.'91 
3.'26/0'.'37 


Bi = 


J=14.41 H=13.78 K s =13.47 
B=19.00 R=17.00 
=18.57 Ri=17.12 B 2 =18.66 R 2 =16.69 1= 


15.91 


IGR J22014+6034 


1 2MASS J22002116+6033420 
USNO-A2.0 1500-08655368 
USNO-B1.0 1505-0321733 


22 h 00 m 21?17 
22 h 00 m 21!17 
22 h 00 m 2r:27 


+60°33'42('06 
+60°33'42'.'ll 
+60 o 33'42'.'44 


8.'09/0 / .'34 

8f09/tf'38 

8.'08/l'.'13 


Bi = 


J=11.66 H=10.99 K s =10.82 
B=15.50 R=13.7 
=15.41 R!=13.61 B 2 =15.19 R 2 =13.38 1= 


12.29 


2 2MASS J22010934+6034088 
USNO-A2.0 1500-08672098 
USNO-B1.0 1505-0322270 


22 h 01 m 09?34 
22' l 01 m 09?29 
22 fe 01 m 09?28 


+60°34'08^'81 
+60°34'08^'12 
+60°34'08'/00 


2.'17/0'.'11 
2.'18/0'.'67 
2.'18/0'.'82 


Bi= 


J=10.62 H=10.39 K s =10.32 
B=12.70 R=11.00 
=12.57 Ri=11.00 B 2 = 11.72 R 2 =10.86 1= 


=10.69 


3 2MASS J22015416+6038094 
USNO-B1.0 1506-0321404 


22 h 01 m 54?16 
22 fe 01 m 53?99 


+60°38W48 
+60°38'09'.'60 


5.'33/0 / .'23 
5?32/l'.'02 


B 


J=11.50 H=ll.ll K s =10.94 
1=14.59 Ri=11.77 B 2 =- R 2 = 12.46 1=11.74 


4 2MASS J22020837+6030425 
USNO-A2.0 1500-08693376 
USNO-B1.0 1505-0322924 


22 h 02 m 08?37 
22 h 02 m 08?28 
22 h 02 m 08?38 


+60°30'42!'50 
+60°30'42;'83 
+60°30'42<'77 


6.'06/0 / .'18 
6.'05/0'.'68 
6.'06/0'.'43 


Bi = 


J=12.59 H=11.95 K s = 11.71 
B=17.50 R=15.50 
=17.37 Ri=15.27 B 2 = 17.19 R 2 =15.38 1= 


=13.54 



a The catalogs are the Two Micron All Sky Survey (2MASS, Cutri et al. 2003) and the United States Naval Observatory (USNO-B1.0 
and USNO-A2.0, Monet 1998, 2003). The astrometric accuracy of 2MASS and USNO catalogs are <0 / /l (Skrutskie et al. 2006) and Of/2 
respectively (Deutsch 1999; Assafin et al. 2001; Monet 2003). 

b The angular distance from the center of the INTEGRAL error circle/The angular distance from the Chandra position reported by T12 



Investigating the Optical Counterpart Candidates of INTEGRAL Sources 



TABLE 4 

Ew Measurements Of Candidate Counterparts To Igr J04069+5042 and Igr J06552-1146 



ICR J04069+5042 



IGR J06552-1146 



SrII 

H<5 

Cal 

H/3 

Cal triplet 



Ha 



A4077 

A4226 

A6103 
A6122 
A6162 



J04064392+5044469 



0.53±0.06 
1.00±0.04 
1.95±0.30 
-1.07±0.04 



J04064872+5039316 

0.65±0.11 
2.42±0.18 
0.94±0.13 
2.8±0.23 

<0.07 

<0.05 

<0.3 
2.41±0.08 



J06545833-1149119 



0.56±0.04 
0.98±0.14 
1.42±0.20 
-0.45±0.05 



Spectral Class 



K5-7 Ve 



F5-8 III-V 



K5-8 Ve 



Note. — EW measurements are given in A units and by convention, positive values denote the 
absorption lines. 



TABLE 5 

Ew Measurements Of Counterparts To IGR J21188+4901 and IGR J22014+6034 



IGR J21188+4901 



J21182288+4906259 J21183906+4858049 



Ha -1.9±0.24 <0.03 



Spectral Class M3-5 Ve K7-M4 V 



IGR J22014+6034 

J22002116+6033420 J22010934+6034088 J22015416+6038094 J22020837+6030425 



Call H & K 


A3934 
A3968 


-1.02±0.51 
-1.36±0.07 






H/3 






4.22±0.01 


* 


Eel 


A5269 




1.07±0.01 




Ha 




-2.15±0.03 


3.60±0.01 


2.26±0.01 


Call triplet 


A8498 




1.82±0.05 


1.25±0.04 




A8542 




2.49±0.03 


2.42±0.02 




A8662 




2.01±0.01 


2.36±0.01 



Spectral Class K7-M4 Ve F5-G5 III-V G III-IV M2-5 Ve 



Note. — EW measurements are given in A units and by convention, positive values denote the absorption lines. 
Since H/3 absorption of J22015416+6038094 is highly blended, its EW measurement is not given in the table. 
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TABLE 6 
Summary Table 



Candidate #Chandra Name 


Spectral Class 


IGR J04069+5042 

Type 


Notes 


iCXOU J040643.9+504446 
2 CXOU J040648.7+503931 


K5-7 Ve 
F5-8 III-V 


Active star; RS CVn? 
Late-type star 


Ha emission, lower SAP a 
No emission lines, higher SAP 


IGR J06552-1146 


iCXOU J065458.3-114911 
2 CXOU J065523. 6-114601 
3 CXOU J065529. 5-114900 


K5-8 Ve 


Active Star; RS Cvn? 
? 
? 


Ho emission, higher SAP 
No optical spectrum, intermediate SAP 
No optical spectrum, low SAP, positive hardness 


IGR J21188+4901 


1 CXOU J211822.9+490627 
2 CXOU J211839.1+485806 


M3-5 Ve 
K7-M4 V 


Active Star; Flare star? 
Late-type star 


Ho emission, higher SAP 
lower SAP but not detected in 2-10 keV 


IGR J22014+6034 


iCXOU J220021. 1+603342 
2 CXOU J220109.3+603409 
3 CXOU J220154.1+603809 
4 CXOU J220208.4+603042 


K7-M4 Ve 
F5-G5 III-V 
G III-IV 

M2-5 Ve 


Active star 
Late-type star 
Late-type star 
Active Star; Flare star? 


Ha emission, high SAP, detected in 2-10 keV 
No features, high SAP 
No features, lower SAP 
Ho emission, higher SAP 



a SAP: Spurious Association Probability 
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Fig. 1 — DSS-R images of the fields around 1GR J04069+5042 (upper left panel), IGR J06552-1146 (upper right panel), IGR J21188+4901 
(lower left panel) and IGR J22014+6034 (lower right panel). Blue circles indicate the 90% confidence radii of INTEGRAL error circle, 
whereas red circles show the positions of the candidate Chandra counterparts in these fields given in T12. For each frame the observed 
candidates which are presented in this work are labeled with numbers. 
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Fig. 2. — The spectrum of candidate #1 of IGR J04069+5042. A weak Ha emission line can be distinguished from the noise level. This 
source is identified as a K5-7 type main sequence star due to the existence of molecular MgH bands and the similarities between the stars 
of the same spectral type. 
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Fig. 3. — The spectrum of candidate #2 of IGR J04069+5042 shows the characteristics of an F type star. Blue arrows indicate the Fel 
lines at AA4046, 4271, 4383, 4819, 4957, 5216 and 5270 respectively whereas symbol ffi shows the telluric absorption band. The Balmcr 
series lines and the CH-G band are the most prominent features. Neutral and singly ionized Magnesium lines, Cal triplet, NaD absorption 
and the blend of several metal lines at 6497 A are also visible. 
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Fig. 4. — R,-band images of the field around the candidates #2 (left panel) and #3 (right panel) of IGR J06552— 1146. The black circles 
show the positions of Chandra sources, but the radii are larger than the error radii for clarity. The Chandra error radii are 0'.'68 and O'.'Sl 
respectively. For the candidate #2, the source in the error circle is ambiguous while the Chandra position is not coincident with any optical 
or IR counterpart in the catalogs for the latter. The i?-band magnitude of the weak source in the error circle for the candidate #3 is found 
to be 22.17±0.14 from the aperture photometry. 




Fig. 5.— The flux-calibrated spectrum of candidate #1 to IGR J06552-1146. Fel lines (red bars) at AA4046, 4271, 4325, 4383, 4404, 
4528, 4871, 4889, 4918, 4957 are shown in a separate panel for clarity. A weak Ho emission is also seen in the main panel. 
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Fig. 6. — The identified features in flux-calibrated spectrum of candidate #1 to IGR J21188+4901. The spectrum contains both broad 
molecular TiO bands with a red continuum and a prominent He* emission, indicating an M type chromospherically active star. The telluric 
absorption bands are denoted by ffi symbol. 
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Fig. 7. — The flux-calibrated spectrum of candidate #2 to IGR J21188+4901. The source has typical features of a late-type star. The 
symbol ® indicates the telluric absorption bands. 
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Fig. 8. — The flux-calibrated spectrum of candidate #1 to IGR J22014+6034. The emission lines of Call H & K and Ha can be considered 
as indicators of a chromospherically active star. 
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Fig. 9.— The flux-calibrated spectra of 2MASS J22010934+6034088, candidate #2 to IGR J22014+6034. The spectrum is dominated 
by Balmer series lines and metallic lines of Fel, Mgl. The symbol © denotes the atmospheric telluric bands. 




Fig. 10.— The flux-calibrated spectrum of 2MASS J22015416+6038094, candidate #3 to IGR J22014+6034. The main spectral features 
Mgl, Mgll and Fel lines in addition to Balmer series lines are labeled. The symbol © denotes the atmospheric telluric bands. This candidate 
is identified as a G type star. 
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Fig. 11. — The flux-calibrated spectrum of candidate #4 to IGR J22014+6034. The spectrum shows typical spectral features of a M 
type main sequence star dominated by strong molecular bands of TiO. 
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Fig. 12. — H — Ks vs J — H diagram for all sources discussed in this work. The validated symbiotics are taken from Phillips (2007), and 
shown with small, empty, purple circles. The solid lines enclose the S type symbiotics, whereas the single dashed line and the first solid line 
enclose D type symbiotics (Corradi et al. 2008). The locus of points for main sequence stars and red giant branch stars are represented by a 
blue solid line and a red dashed line respectively. The questionable symbiotic counterparts (filled brown symbols) are for IGR J16358— 4726 
and IGR J17497-2821, the incorrect IGR symbiotic candidates (empty red circles) are for IGR J11098-6457 and IGR J17197-3010, and 
the incorrect IGR symbiotic identified through K-haxi& spectroscopy is for IGR J16393— 4643. See §4.5 for more details. 



